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A superfluid-insulator transition is known to occur in strongly disordered Fermi gases, in both 
the BCS and BEC regimes; here, we address the properties of this transition across the BEC- 
BCS crossover. We argue that the critical disorder strength at which superfluidity is lost changes 
non-monotonically with detuning from Feshbach resonance, and that a reentrant superfluid phase 
arises for detunings near the fermionic mobility edge. Our analysis of the intermediate regime is 
quantitatively vahd for narrow resonances and near four dimensions, and provides a simple physical 
picture of this regime, in terms of two distinct but coexisting insulators. 
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I. INTRODUCTION 

The physics of interacting electrons in strongly disor- 
dered media involves a delicate interplay between Cooper 
pairing, Anderson localization, and Coulomb repulsion. 
The nature of the resulting low-temperature phases has 
long been an active topic of research [IHZj: although 
a quantum superconductor-insulator transition (SIT) is 
generally believed to occur, the actuating mechanism(s), 
the nature of the insulating phases and of their charge 
carriers, and the possible existence of a zero-temperature 
metallic regime [21 E] are still open and much-discussed 
questions [7H5]. Theories of the transition vary, e.g., in 
whether they incorporate the composite nature of Cooper 
pairs and the repulsive and/or long-range parts of the 
Coulomb interaction, and in the nature of the disorder 
considered; the importance of each of these factors in 
the condensed- matter systems studied is not, in general, 
understood a priori. 

It is thus hoped that the microscopically well- 
characterized settings provided by ultracold atomic 
gases, trapped in controllably disordered environ- 
ments pm [TTj. can shed light on the SIT. Beyond this 
general motivation for studying the SIT in ultracold 
atomic gases, we note that two major approaches to the 
SIT correspond to the limits of BCS and BEC pairing 
respectively — namely, (i) those that investigate the BCS 
pairing of localized electrons, ignoring phase fluctuations 
and repulsion [TJ |8] ; and (ii) those that address tightly 
bound, repulsively interacting bosons via (e.g.) the Bose- 
Hubbard model, ignoring their composite character [3]. 
As these approaches correspond respectively to the BCS 
and BEC limits, interpolating between them is equivalent 
to understanding how the SIT evolves as one tunes across 
the BEC-BCS crossover via a Feshbach resonance [T^ . 
This question has not, to our knowledge, been previ- 
ously considered in detail; the limiting cases are discussed 
in Ref. [T3], and disorder effects on the superconduct- 
ing transition temperature have been treated perturba- 
tively [14 . As these works observe, superconductivity is 
more robust to disorder in the BCS regime than in the 
BEC regime, owing to Anderson's theorem |T5]; how- 



ever, the properties of the SIT in the intermediate, near- 
unitary regime have received little attention. 

In the present work we address the nature of the 
SIT in the intermediate region, focusing on the mor- 
phology of the zero-temperature phase diagram. Our 
treatment is quantitatively valid for narrow Feshbach 
resonances and/or dense gases [16l [17]; the former of 
these limits can be realized in experiments with ^Li [TS] 
(which has a narrow resonance at approximately 543 
G [Tni[in]). In the narrow-resonance limit, we show that 
the superconductor-insulator phase boundary is, for a 
certain range of disorder strengths, nonmonotonic as a 
function of detuning from the Feshbach resonance, in- 
volving an intermediate insulating phase for detunings 
slightly below the mobility edge. We establish this cen- 
tral, unexpected, result both for bounded disorder and 
for unbounded Gaussian disorder (which is most relevant, 
e.g., to optical speckle potentials). We argue that the 
shape of the phase diagram can be understood in terms 
of a competition between the pairing instability and the 
tendency of Cooper pairs, which are bosons, to condense 
into deeply localized states; thus, the intermediate in- 
sulating phase has fermions of large localization length 
coexisting with deeply localized bosons. We discuss the 
feasibility of realizing this intermediate insulating phase, 
and finally note possible connections between our anal- 
ysis and the e-expansion approach pioneered by Nishida 
and Son [IDj. 



II. MODEL AND APPROACH 



We consider the two-channel model of a Feshbach res- 
onance |16) . This model involves fermionic gases con- 
sisting of consisting of two "spin" species (i.e., internal 
states), coupled to a gas of bosonic molecular states, such 
that two fermions can pair up to form a bosonic molecule: 
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This Haniiltonian was analyzed for clean systems in 
Ref. |16j in the small-g limit; note, however, that a 
weakly-coupled two-channel description applies near uni- 
tarity for all short-range fermionic interactions near four 
dimensions [20) . For disordered systems, we have sup- 
plemented the Hamiltonian of Ref. [16] with a contact 
repulsion of strength U among the bosons, which is as- 
sumed to be smaller than the other energy scales; this 
contact term, always physically present, can be neglected 
for clean systems, but must be included for disordered 
systems in order to circumvent the pathological behavior 
of the disordered noninteracting Bose gas. As we shall 
see, our results are insensitive to U. The potential T^(x) 
includes both disorder and (in principle) a periodic po- 
tential (note, however, that the narrow-resonance model 
does not correspond to the attractive-C/ Hubbard model, 
as pair hopping is explicitly included in H and does not 
arise by a superexchange mechanism). We shall focus 
on three- (or higher-) dimensional systems, in which the 
single-particle (i.e., 5 = 0) Hamiltonian has both local- 
ized and extended eigenstates, separated by a mobility 
edge [H] at some energy Ef.- For most of the present 
work, we shall take the disorder distribution to be speci- 
fied through E(, (which scales with the disorder strength) 
and the single-particle density of states p{E). We first 
analyze the case of a bounded disorder distribution with- 
out spatial correlations (e.g., a box distribution [21]) and 
then turn to the case of unbounded Gaussian disorder 
with short-range correlations (which is relevant, e.g., to 
the current speckle-potential-based experimental realiza- 
tions [llj). While our work does not explicitly treat the 
case of smoothly varying or correlated disorder, such an 
extension would be straightforward. 

Our approach is as follows. First, we consider 
the Hamiltonian H in the limit g = Q (known as 
the infinitely-narrow resonance limit [16j). in which it 
consists of decoupled fermionic and bosonic channels 
(Fig. [T]). The equation of state for each of these chan- 
nels is known, and these results can be combined to form 
a "zeroth-order" theory valid throughout the crossover. 
Having established the structure of this 3 = theory, we 
include the effects of g perturbatively, to lowest order, 
adapting the results of Ref. [5] . This perturbative analy- 
sis shows the existence of regimes of reentrant supercon- 
ductivity, which, we argue, offer a conservative estimate 
of the parameter range in which reentrant behavior takes 
place. The primary motivation for choosing the narrow 
resonance approach is that it enables us to study the en- 
tire BEC-BCS crossover while working at strong disorder 
(in contrast with previous works [14j). 





FIG. 1 : Limit of a narrow Feshbach resonance, (a) Deep in the 
BCS limit, the bosonic channel is above Ef and is therefore 
unoccupied; all the particles are in the Fermi channel. The 
Fermi energy Ef is taken to be above the mobility edge Ec- 
(b) Crossover region: e is lower than Ef\ therefore, fermions 
with energies above e pair and migrate into the molecular 
channel. The molecules cannot all populate the lowest bosonic 
state owing to repulsive interactions between them; therefore, 
they occupy a finite energy range. 



III. LIMITING CASES 

In this section we briefly outline the general properties 
that the SIT is expected to exhibit in the BEG and BGS 
limits, without specific reference to the narrow-resonance 
model. For the purposes of this section we shall assume 
that there is no underlying lattice. We then proceed, in 
the following sections, to show that the narrow-resonance 
approach both reproduces these expectations and pro- 
vides a useful scheme for interpolating between them. 

BCS limit. In the BGS limit, the nature of the ground 
state depends on the location of the Fermi energy Ep 
with respect to the mobility edge i?c [2 HI US] • In par- 
ticular, if Ep > Ec then for any infinitesimal attraction 
(i.e., any detuning) the ground state is superconduct- 
ing 1_. On the other hand, if Ep < Ec then a finite 
minimum attraction [8] is necessary to give rise to a su- 
perconducting state. Thus, for Ep < Ec an insulating 
ground state obtains deep in the BCS regime, but as the 
detuning is decreased (and the BGS coupling becomes 
stronger) a transition to a superconductor takes place. 
Thus, at fixed disorder, the critical density behaves as in 
the "BGS" side of Fig. [Sjja): it is asymptotically given 
by the criterion Ep = Ec, but closer to resonance the 
critical density decreases. 

BEC limit. In the BEG limit, the SIT is governed by 
a competition between the repulsion between the Gooper 
pairs and the disorder strength. For weak repulsive inter- 
actions the bosons occupy isolated deep wells of the dis- 
order potential; as the repulsive interaction strength (or 
density) is increased, they occupy more and more wells 
of the disorder potential, until eventually the Josephson 
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coupling between the wells is strong enough to give rise to 
global phase coherence. An approximate criterion [T3t[22] 
for the critical repulsive interaction strength is given by 
Ec « nab/{4:M), where Ub is the boson-boson scattering 
length (which decreases to zero as one tunes away from 
resonance). Thus, sufficiently deep in the BEC limit, the 
ground state is always insulating, as shown in the "BEC" 
side of Fig. |3ja); this is in sharp contrast with the BCS 
limit discussed above. 

Crossover. The considerations above suggest that, for 
a given disorder strength, the lowest critical density for 
the SIT occurs in the crossover regime (or, equivalently, 
that superconductivity is most robust against disorder in 
the crossover regime). Thus, the critical disorder is ev- 
idently a non-monotonic function of density, with quite 
different functional forms in the two limits. In what fol- 
lows, we shall discuss this nonmonotonic function at a 
more quantitative level. 



IV. CASE OF BOUNDED DISORDER 
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FIG. 2: Change of the chemical potential /i with detuning, 
in the infinitely narrow resonance limit for finite U, in the 
case of bounded disorder (red dashed Une) and unbounded 
Gaussian disorder (black solid line). Energies are expressed 
relative to the Fermi energy. The two cases are similar in the 
BEC regime but differ in the BCS regime. 



We first address the conceptually simple case of 
bounded disorder, and then extend our results to un- 
correlated Gaussian disorder. As discussed above, we 
consider the decoupled {g = 0) limit first, then work per- 
turbatively in g. 



trast with the clean case [TB], we must self-consistently 
determine the chemical potential, as the bosons trapped 
in deep rare wells interact. We first determine the chem- 
ical potential due to a dilute gas of interacting bosons at 
density n^, adapting the analysis of Ref. [32] as detailed 
in the Appendix; the result is: 



A. Infinitely narrow resonance limit 

We begin with the limit of an infinitely narrow res- 
onance (i.e., g —7- 0), so that the fermions and bosons 
are decoupled. (Following Ref. [TB], we assume that the 
g — > limit is taken in such a way that the two chan- 
nels equilibrate with each other.) For e 3> Ep (i.e., in 
the BCS limit), the Fermi sea is filled up to a Fermi 
energy Ep; depending on whether this exceeds E^ the 
system is either an insulator or a metal. On the other 
hand, for — e ^ Ep (i.e., in the extreme BEC limit), 
all fermions are tightly bound into molecules interacting 
via the weak contact repulsion U. The case of U = 
is pathological, as it involves all bosons occupying the 
lowest single-particle state; for small nonzero U, how- 
ever, the bosons form a "Bose glass" consisting of widely 
separated, mutually phase- incoherent, puddles of bosons 
at deep disorder minima: thus, the system is insulat- 
ing regardless of the relative locations of Ep and Ec- 
We shall assume that U (or the corresponding scattering 
length oo) is small enough, and the relevant densities low 
enough (i.e., nao/2m ^ Ec), that the molecules are lo- 
cated in the exponentially rare Lifshitz tail states "21 of 
the disorder potential. 

In the g — J> limit, the crossover between these limits 
occurs as follows. As one decreases e until it is lower than 
Ep, the chemical potential fi begins to track e, as it is fa- 
vorable for fermions with energies greater than e to pair 
up into molecules, which then form a Bose glass. In con- 



In 



no 
rib 



-1 -2/3 



(2) 



where £ is the characteristic energy scale associated with 
the disorder, and ng ^ where ab is the scatter- 

ing length for bosons. (Note that the chemical potential 
therefore depends weakly on U.) Given this expression, 
one can determine the density of bosons in the interme- 
diate regime, assuming they are able to equilibrate with 
the fermions, by setting the two chemical potentials to be 
equal, and noting that the number of molecules created is 
one half the number of fermions depleted from the Fermi 
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where p{E) is the density of states for the fermions. Note 
that Eqs. ([2] [3]) together comprise an equation of state 
for the composite system in the g = limit; the depen- 
dence of the chemical potential on detuning is shown in 
Fig. [2j The qualitative behavior they imply is as fol- 
lows. In the "metallic" case where Ep > Ec, n decreases 
past Ec for some e; at this point, the fermionic sector 
undergoes a metal-insulator transition, and consequently 
(as the bosons are trapped in isolated wells) the com- 
posite system becomes insulating. On the other hand, in 
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the "insulating" case where Ep < Ec, the system is al- 
ways insulating. Therefore, in this limit one is led to the 
schematic phase diagram denoted by the black dashed 
lines in Fig. [sf^a). 

The above discussion extends trivially to the case of a 
lattice with one or more bands; each band generically 
has an upper and a lower mobility edge, and metal- 
insulator transitions occur whenever the chemical poten- 
tial (which is approximately e) coincides with a mobility 
edge. For sufficiently small [/, bosons always occupy lo- 
calized states near the bottom of the lowest band, and 
do not contribute to transport. 

B. Finite-width resonances: BCS and BEC limits 

We now consider how the above picture is altered for 
small but finite g. We first briefly sketch how the narrow- 
resonance formalism reproduces the known qualitative 
features of both the BCS and BEC regimes [T3], and 
then turn to the intermediate regime. 

a. BCS regime, e 3> Ep. In this regime, there is a 
gap to the creation of molecules; therefore, the bosonic 
channel can be integrated out. The nature of the ground 
state depends on the position of Ep relative to Ec- For 
Ec < Ep, the Fermi energy is among the extended states 
(i.e., the g = state is metallic), and BCS pairing occurs 
for arbitrary g > and arbitrary e, with a gap scale set 
by A ~ Epexp[-e/{p{Ep)g^)]. For Ec > Ep, the Fermi 
energy is among the localized states; pairing occurs in 
this regime, but only for couplings above some critical 
value gc (or, more relevantly to the present discussion, 
above some critical value Sc), for which a rough estimate 
is given by the following implicit expression [8] |9| : 

i^n^) ^-p[^^/(^(^^)^'')] - 1- (4) 

where i' is the localization length critical exponent, which 
is approximately unity |21j . Thus, at zero tempera- 
ture the system is insulating in the deep BCS regime 
(e — > +oo), but undergoes an SIT upon sweeping toward 
resonance, as s decreases and the effective superconduct- 
ing coupling g'^ /e correspondingly increases. 

b. BEC regime, e ^ 0. In this regime, all fermions are 
dimerized and there is a gap to the breaking of molecules; 
therefore, the fermions can be integrated out, giving rise, 
for any g > 0, to an effective repulsive interaction among 
the molecules. The strength of this induced interaction 
depends on the associated energy denominator. In par- 
ticular, deep in the BEC regime, the induced interaction 
is parameterized by the scattering length |16j 

which decreases as e — >■ — oo. Note that, while the 
narrow-resonance behavior of ai, differs quantitatively 



from the broad-resonance value [TB], the trend that in- 
teractions weaken far from the resonance is true in both 
cases. Moreover, as previously discussed, our results for 
the disordered system depend only logarithmically on ab, 
so that the distinction between the narrow- and broad- 
resonance scattering lengths is immaterial for our pur- 
poses. Note that if the quantity 

4m 

exceeds the characteristic disorder scale £, the Bose glass 
undergoes a transition to a supcrfluid state. Conse- 
quently, the ground state in the BEC regime is a su- 
perfiuid for sufficiently high densities and interaction 
strengths, whereas in the BCS regime it is a superfluid 
at high densities, regardless of the interaction strength. 

C. Intermediate regime, < e < Ep 

The intermediate regime poses greater difficulties: in 
this regime, the system possesses both fcrmionic and 
bosonic excitations at the chemical potential, so that, in 
principle, neither channel can safely be integrated out. 
The argument of the present section is that one can nev- 
ertheless compute certain properties of the phase dia- 
gram in this regime, and show, in particular, that the 
critical coupling is a nonmonotonic function of e. For 
specificity, we shall focus during this argument on the 
case of a three-dimensional Anderson lattice model [21] , 
with on-site energies drawn from the uniform distribution 
[—Wt/2, Wt/2], where t is the nearest-neighbor fermionic 
hopping matrix element; however, the results are evi- 
dently more generally applicable. (The bosonic hopping 
is ^ t/2 for this model; note, again, that the two-channel 
model on a lattice is distinct from the attractive-[/ Hub- 
bard model.) 

We assume that W is close to but not quite at its criti- 
cal value [21], Wc ~ 16.5, so that Ec is near the center of 
the band; we assume, further, that Ep is slightly above 
Ec- The density of states p{E) is known [23] to be ana- 
lytic around Ec, so it follows that, if Ec is near the center 
of the band, p{E) is approximately constant near Ec- Un- 
der these conditions, the ground state is superfiuid for in- 
finitesimal g in the BCS regime, as previously discussed; 
however, as e is decreased until fi < Ec, the ground state 
becomes insulating (in both the fermionic and bosonic 
sectors) in the decoupled limit, 5 — 0. For g 0, there 
are two kinds of low-energy states: (a) fractal, almost- 
extended fermionic states, and (b) deeply localized Lif- 
shitz tail [21] [25] states, each occupying a rare fluctua- 
tion of the disorder potential and holding a small number 
of bosons (these become arbitrarily rare as Ep — ?> Ec). 
As both channels are stable against superfluidity for in- 
finitesimal g, a minimum inter-channel coupling gc > 
is required for global superfluidity to occur; and under 
the assumptions above it is clear that, at T = 0, global 
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superfluidity first arises via a fermionic pairing mecha- 
nism as g is increased through gc- We argue that, under 
the assumptions above, gc is in general a non-monotonic 
function of fi and therefore of s. 

We now discuss how pairing arises in this regime, 
and estimate the pairing strength. In contrast to the 
BCS and BEC limits, the intermediate regime requires 
us to consider which bosonic states mediate the pair- 
ing interactions — in the clean system, one can assume 
that pairing is mediated by the zero-momentum bosonic 
state , but this has no direct analog in the disordered 
system. To address this issue, we first note that there 
are, in principle, two kinds of bosonic states that can 
mediate pairing: (i) the exponentially rare, but strongly 
interacting Lifshitz tail states within g of the chemical po- 
tential, and (ii) the more abundant but relatively weakly 
interacting states with appreciable gaps. We argue that, 
under the assumptions listed above, the Lifshitz-tail con- 
tribution to pairing is negligible, so that all the relevant 
bosonic degrees of freedom are in fact gapped and can 
be integrated out, allowing us to adapt the BCS-limit re- 
sults with some effective energy denominator Ej^t- (We 
shall qualitatively address the role of the tail states in 
the next subsection, and show that they accentuate the 
nonmonotonicity of the phase diagram.) 

The argument for neglecting the Lifshitz tail states 
runs as follows. The resonant pairing mediated by the 
tail states within a window dE of the chemical potential 
can be estimated as 

M,:cs. ~ g^/nbPb X overlap (7) 

where Uh is the density of bosons; Pt, is the number of 
resonant tail states ~ exp[— — e))^'^/^], where B is 
related to f^^^; and the coupling is modulated by a typi- 
cal overlap integral (of the form J IV'P^) associated with 
the fermion-boson coupling [Eq. ([T])]. By contrast, the 
off-resonant pairing mediated by higher-energy bosonic 
states within an energy window dE can roughly be esti- 
mated as 

M{E) — Pb(E)dE X (overlap)^ (8) 

E — fi 

where pi, is the bosonic density of states at E. Compar- 
ing these two expressions, one finds that (even ignoring 
the nature of the overlaps), the exponential increase in 
the number of states outweighs the decreased coupling 
to each state, so long as g/{Wt) >• y/NlPij{fJ.) , which is 
always true for fi w Ep (i.e., Ni, — > 0). Furthermore, 
the numerical results of Ref. [Mj suggest that wavefunc- 
tions far apart in energy are in fact spatially anticor- 
related. This effect would further suppress the interac- 
tions between the fermions and the tail-state bosons, and 
therefore extend the regime of quantitative validity of the 
present analysis. 

Having established that the tail states contribute neg- 
ligibly to pairing, we now estimate which states do pro- 



vide the dominant pairing interactions. One can eas- 
ily see that states that are too high in energy cannot 
contribute to pairing, as for these states the matrix ele- 
ments of the form J dx\tp{x)\'^ip{x) are products of three 
rapidly oscillating functions and must vanish. Combin- 
ing these observations (viz. that both states with very 
small energy denominators and those with large energy 
denominators are negligible), one sees that M{E) must 
be peaked at some E^t measured from s. This en- 
ergy Elt is of order Wt in general; for the case of the 
three-dimensional Anderson model, the considerations 
of Ref. |i24 (regarding the "stratification" of coordinate 
space) would suggest a value Elt ~ s + Wt/2 [3T]. How- 
ever, for our purposes we only need the fact that it is 
nonzero, which enables us to estimate a BCS pairing in- 
teraction of strength 7(/i) ^ g'^ / {E LT — p) ■ It is crucial to 
note that j{fJ.) increases as one tunes toward resonance, 
because -Elt = e+ constant, and fi ~ e increases as the 
tail states fill up with bosons. 

The arguments above lead us to a problem involving 
only fermionic states at the chemical potential and vir- 
tual bosonic states ~ Ej^t away from it. This reduced 
problem can now be addressed exactly as in the BCS 
limit, by inserting 7(/Lt) into Eq. Q, to arrive at the fol- 
lowing criterion for superfluidity: 



[~\^) /^(M)7(M)exp[l/(p(A.)7)]<l. (9) 

In this expression, the relation between p and e is given 
by Eqs. ([2)|3|. Thus, Eq. ^ thus deter mines a critical 
value of 7 and thus of g, which is plotted in Fig.jsjjd), and 
can be seen to be nonmonotonic, exhibiting a transition 
to an insulator very close to E^. as well as a reentrant 
superfluid at lower detunings [Fig. [3][c)-(e)]. The intu- 
itive reason for this non-monotonic behavior is that the 
bosonic tail states fill up rapidly (there being very few 
such states), and hence decrease the energy denomina- 
tor for pairing (Elt — p) without substantially depleting 
the Fermi sea. Note that as e is swept closer still to 
resonance, the fermionic density of states at p rapidly 
decreases, triggering a second transition to an insulating 
state, which persists into the BEC limit. This second 
transition generally occurs far from the mobility edge, 
and is thus beyond the scope of the present work. (How- 
ever, it follows from our considerations and the previ- 
ously mentioned results for the BEC limit |13j that such 
a transition must occur at some point.) 

Let us briefly review the argument given above for 
the shape of the phase diagram. Consider a situation 
in which Ep is slightly above Ec, and the detuning e is 
tuned so that p is slightly below Ec- The ground state 
is an insulator for g = 0, and this insulator is evidently 
stable for sufficiently small, nonzero g; however, it be- 
comes a superconductor at some gc- This critical cou- 
pling gc depends on p, (or equivalently on e) in three 
ways: (1) the localization length at p decreases as p de- 
creases, suppressing pairing; (2) the fermionic density of 
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states decreases as ^ decreases, suppressing pairing; 
and (3) the pairing interaction grows stronger as /i de- 
creases, owing to the concomitant decrease in the energy 
denominator Et — /i- Sufficiently close to the mobility 
edge, effect (1) dominates; sufficiently deep in the insu- 
lating phase, effect (2) dominates; in either case pairing 
is suppressed with decreasing /x. However, as argued in 
this work, there is a range of intermediate energies for 
which effect (3) is the dominant effect, so that pairing 
is enhanced by a decrease of leading to a nonmono- 
tonic phase diagram overall. (The competition between 
effects (2) and (3) has previously been noted in the 
condensed-matter context, as leading to a maximum in 
the transition temperature of SrTiOa as a function of 
carrier density.) 



D. Role of the Lifshitz tail states in mediating 
pairing 

In the preceding discussion, we neglected the role 
played by the Lifshitz tail states in mediating pairing, ar- 
guing that this role was negligible for Ep w Ec- We now 
return to this point. The role of the Lifshitz tail states is 
to provide an additional channel for pairing, of strength 
~ y/nl where rif, is the density of molecules. Now, as e is 
decreased, the bosonic band fills up, so that both and 
rif, increase. Thus, the coupling strength increases faster^ 
as £ (or equivalently /i) is reduced, than our previous 
analysis predicted. This strengthens effect (3) discussed 
in the previous paragraph, and thus accentuates the non- 
monotonicity of the phase diagram. 



E. Phase diagram 

The considerations in the previous sections can be syn- 
thesized into a global phase diagram, at fixed disorder 
strength, in terms of the density, detuning, and coupling 
g (i.e., the resonance width). The results are sketched 
in Fig. [3]Ja)-(c). In the narrow-resonance limit, there 
are four possible scenarios for the BEC-BCS crossover, 
as follows. (1) When Ep <C E^ the system stays in- 
sulating through the crossover. (2) When Ep is near, 
but lower than E^, superfluidity occurs at some finite, 
positive detuning e, before vanishing at yet smaller e. 
(3) When Ep is slightly higher than i?c, superfluidity is 
destroyed for e just below the mobility edge, but reap- 
pears as the detuning is decreased further, before finally 
vanishing again in the BEC limit. (4) For Ep ^ Ec, 
the ground state is superfluid in the BCS and intermedi- 
ate regimes, and superfluidity persists well into the BEC 
regime before eventually disappearing deep in the BEC 
regime, as the interactions between bosons weaken. Of 
these scenarios, the most exotic is Scenario (3), and one 
might ask to what extent it is limited to the narrow- 
resonance approach adopted here. The crucial point in 
our arguments above where the narrow-resonance limit 
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FIG. 3: (a) Schematic T = phase diagram of the SIT as 
a function of detuning. Main panel shows the phase dia- 
gram at fixed disorder strength as density (i.e., Ep) is varied: 
black dashed line corresponds to an infinitely narrow reso- 
nance, whereas the brown and red lines correspond to rela- 
tively narrow and broad resonances respectively. is the 
mobility edge. The results of this paper concern the shaded 
region. Inset shows phase diagram as a function of disorder at 
fixed density; color coding is the same as in the main panel. 

(b) Sketch of Tc as a function of detuning near an SIT; the 
proximity to an insulating phase is expected to depress Tc. 

(c) -(e) Inverse of the left-hand side of Eq. ([9| as a function 
of ^. Panel (c) shows the case of bounded disorder at three 
densities, from Ep = 1.18i?c (top) to Ep = l.lEc (bottom), 
at fixed r = p{Ec)g^ / Ec = 1/20; panel (d) shows the case 
of Gaussian disorder, with p{Ec)g^ / Ec — 1/3; and panel (e) 
shows the case of Gaussian disorder, with Ep = l.lEc and 
different values of P. The insulating region is shaded gray. 



was invoked was in justifying non-degenerate perturba- 
tion theory for the pairing interaction in the intermediate 
regime; while strictly valid only for g <^ Epp ~ E^, one 
expects that it will remain reasonably accurate so long 
as (7 < Ec- (By contrast, in the more commonly studied 
broad-resonance limit, g 3> Ep « Ec, so that the varia- 
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tion of the energy denominator is irrelevant and one does 
not expect a re-entrant superconducting phase.) How- 
ever, this limitation can be circumvented either by work- 
ing at high densities and strong disorder, or by choosing 
a narrow Feshbach resonance (as discussed below); be- 
sides, even beyond its regime of quantitative vahdity, the 
underlying mechanism by which strong pairing is offset 
by the tendency of bosons to localize might lead, e.g., 
to anomalous behavior in the critical temperature of a 
disordered Fermi gas near unitarity. 



F. Nature of insulator 

The discussion above suggests a simple physical picture 
of the insulating phase in the intermediate regime: viz. 
that it consists of two interpenetrating insulators, a Lif- 
shitz glass (TUl [22] and a "fractal insulator" [S] , which 
are essentially mutually noninteracting because of the 
lack of wavefunction overlap. Thus, the gap and other 
properties of the insulator — measured, e.g., via rf spec- 
troscopy j27j — near the SIT should consist of an appro- 
priately weighted average of the properties of these two 
systems. The coupling between the systems grows for 
smaller e and/or higher densities, as the chemical poten- 
tial then no longer lies in the bosonic Lifshitz tail (de- 
fined as the region in which the density of states grows 
exponentially with energy); in this regime, a more suit- 
able model is to treat the fermions as being strongly 
proximity-coupled to puddles of bosonic superfluid [6] . 



V. UNBOUNDED GAUSSIAN DISORDER 

We now extend our results to the case of uncorrelated 
Gaussian disorder of the same characteristic energy scale 
£, the case most directly relevant to recent experiments 
in optical lattices HI]. The essential distinction between 
bounded and unbounded disorder is the nature of the Lif- 
shitz tail [21]; in the uncorrelated Gaussian case, the den- 
sity of states in the tail scales as p{E) ~ exp{—B\E\-^/^). 
For unbounded disorder, there is some finite probabil- 
ity of finding states at arbitrarily low energies, so the 
bosonic band is not bounded below. Thus, the BCS 
and BEC regimes are not sharply separated: at any de- 
tuning, there are some molecular states as well as some 
fermionic states, and the system is always in the inter- 
mediate regime discussed above. To extend the previous 
analysis to this case, it suffices to replace Eq. ^ with 
the following expression |22) : 



^ — e ~ £ 



In 



Ep 



(10) 



Note that Eq. (10 1, unlike Eq. ([2|, predicts that, regard- 
less of e, /i < Ep (although the difference vanishes as 
£" — in the clean system limit). This replacement 



changes various quantitative features of the phase dia- 
gram, but our qualitative conclusions (and in particu- 
lar, the possibility of a nonmonotonic superfiuid-insulator 
boundary) still hold (see Fig. [3f^d)-(e)). 

Another frequently considered, unbounded disorder 
distribution is the Lorentzian distribution, for which the 
density of states can be computed exactly [28^; how- 
ever, this distribution does not possess an exponential 
Lifshitz tail (instead, the tail density of states falls off 
as p{E) ~ l/E"^). Thus, the relative contribution of 
tail states to pairing is much larger than in the bounded 
or Gaussian cases. While the qualitative considerations 
of the previous section are expected to apply here, the 
parameter-dependence of the transition line is expected 
to differ substantially. 



VI. DISCUSSION 

This work has addressed the behavior of the superfiuid- 
insulator transition for disordered Fermi gases across a 
BEC-BCS crossover, focusing on the intermediate (i.e., 
near-unitary) regime. We have established that the 
superfiuid-insulator phase boundary is nonmonotonic as 
a function of detuning, so that an experiment sweeping 
across a Feshbach resonance can in principle see three 
separate transitions. We were able to establish these 
conclusions quantitatively only in the narrow-resonance 
regime (i.e., where the energy scale associated with 
the Feshbach resonance width [T^ is smaller than the 
Fermi energy); however, it is plausible that, as long as 
the resonance is not too broad (i.e., if g ~ Ep), the 
proximity to an insulating phase will manifest itself in 
a nonmonotonic critical temperature. Moreover, nar- 
row Feshbach resonances do exist (e.g., g/Ep 0.1 in 
^Li [HJ [IS]), and have been proposed as a way of real- 
izing higher-temperature superconducting transitions in 
ultracold atomic systems [T^. Finally, one should note 
that, as the width of the resonance is measured relative to 
the Fermi energy, all experiments at sufficiently high den- 
sities (e.g., in tightly confining traps) can be described by 
a narrow-resonance model. In order to stress the impor- 
tance of the mobility edge, we have restricted our anal- 
ysis to three-dimensional disordered systems; however, 
the generalization to systems of reduced dimensionality 
is straightforward. 

A separate limitation of our analysis, which ap- 
plies generally to the dynamics of disordered quantum 
systems, is that in considering the equilibrium state 
we have neglected the experimentally relevant ques- 
tion of whether equilibrium is achieved on the relevant 
timescales (which must, at a minimum, be much larger 
than 1/g in the narrow- resonance limit). Indeed, recent 
work on the many-body localization problem |29j sug- 
gests that closely related systems do not equilibrate in 
the thermodynamic limit. 

Near-unitary system in 4 — e dimensions. Finally, we 
note that the weakly-coupled two-channel model [Eq. ([I])] 
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becomes valid for any Feshbach resonance, near unitar- 
ity, as the dimensionality of space approaches four. (The 
resemblance arises because the two-body wavefunction at 
resonance is sharply peaked when the particle positions 
coincide, so that weakly bound or even resonant fermion 
pairs behave like pointlike bosons [20].) It follows that 
the ground state at unitarity is insulating for arbitrar- 
ily weak disorder near four dimensions. It is straight- 
forward to see that the ground state near unitarity is 
also generically insulating near two dimensions, by the 
following logic. Near two dimensions, unitarity corre- 
sponds to vanishingly attractive interactions; however, a 
finite disorder strength would localize all single particle 
states in two dimensions, thus overcoming an infinitesi- 
mal pairing term. As a consequence, the phase diagram 
and transition temperature of a disordered unitary Fermi 
gas are nonmonotonic as a function of dimensionality, in 
contrast with the clean-system limit [20 . That the SIT 
occurs for weak disorder in both the 4 — e- and 2 + e- 
dimensional limits suggests the e-expansion as a promis- 
ing strategy for future studies of the unitary-gas SIT us- 
ing renormalization-group methods |30j . 
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Appendix: Chemical potential in the bosonic channel 

We briefly sketch how the density of bosons confined 
in Lifshitz tail states can be estimated, for the case of 
bounded disorder. (The unbounded case is addressed in 
Ref. t22j.) For concreteness we assume the on-site disor- 
der takes on two values, ±Eq. Suppose the chemical po- 
tential is in the Lifshitz tail; then, according to a standard 
argument [21] , the single-particle states that can be filled 
must have dimensions of at least i?_ = ■\/2to(/x — Eq). 
In the Thomas-Fermi approximation (valid for the rel- 
atively large wells in the Lifshitz tail), a well of size 
R < R- typically holds Nr w (/x — Eq)R^ /g parti- 
cles; put differently, a typical well at energy E holds 
Ne ^ in- Eo){E ~ Eq)-^/'^ particles. The total number 
of particles is given by 



N = j N{E)p{E)dE, (A.l) 

where we can use the well-known result [H] for the den- 
sity of states, p{E) - exp{-B\E-Eo\~^/^). This integral 
is dominated by the contribution from energies around /i. 
Therefore, to logarithmic accuracy, one can invert the re- 
lation between N and fi to write 



p = Eo + const. [In no /n] (A. 2) 

for some constant, which, on dimensional grounds, must 
be the energy scale £ associated with the disorder 
strength. 
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